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1 .  INTRODUCTION 


Interception  of  hostile  communications  is  attempted  for  many  diverse 
reasons,  such  as  reconnaissance,  surveillance,  position  fixing,  identi¬ 
fication,  or  a  prelude  to  jamming.  Different  purposes  require  different 
systems,  but  whatever  the  purpose,  an  interception  system  nearly  always 
must  achieve  the  three  basic  functions  of  detection,  frequency  estima¬ 
tion,  and  direction  finding.  Although  these  three  elements  of  intercep¬ 
tion  are  usually  integrated  in  a  practical  system,  they  are  discussed 
separately  in  this  paper  for  clarity  of  presentation.  The  basic  con¬ 
cepts  and  issues  of  the  three  elements  are  presented  at  the  systems 
level,  assuming  that  little  is  known  about  the  signals  to  be  inter¬ 
cepted.  Primarily  because  of  the  rapidly  changing  technological  base, 
the  implementation  and  the  engineering  details  of  the  interception 
systems  are  not  addressed.  Although  this  paper  is  concerned  with  the 
interception  of  communications,  only  slight  modifications  of  the  results 
are  required  to  apply  them  to  the  interception  of  radar. 

The  potential  interceptor  has  at  least  one  major  advantage  over  the 
communicators.  The  accuracies  of  detection,  frequency  estimation,  and 
direction  finding  are  determined  by  the  energy  of  the  entire  message 
transmitted,  which  may  include  many  symbols.  In  contrast,  the  intended 
receiver  makes  decisions  with  accuracies  determined  by  the  energy  of 
each  transmitted  symbol.  From  another  point  of  view,  the  intended 
receiver  generally  must  make  many  separate  decisions,  whereas  the  inter¬ 
ception  receiver  must  make  only  a  few  decisions. 


2.  DETECTION 

If  the  form  and  the  parameters  of  the  signal  to  be  intercepted, 
s(t),  were  known,  optimum  detection  in  white  Gaussian  noise,  n(t),  could 
be  accomplished  by  a  matched  filter  or  an  ideal  correlator.  Figure  1 
depicts  a  correlator  for  the  received  signal,  r(t)  =  s(t)  +  n(t),  and  an 
observation  interval ,  T.  The  comparator  input  is  compared  with  a  fixed 
threshold  level,  VT,  to  determine  the  presence  of  an  intercepted 
signal.  It  is  a  standard  result  that  the  probability  of  false  alarm, 
Pp,  and  the  probability  of  detection,  PD,  are  given  by1 


-*A.  Whalen,  Detection  of  Signals  in  Noise,  Academic  Press,  Inc.,  New 
York  (1971). 
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where  E  is  the  signal  energy,  Ng/2  is  the  noise  power  spectral  density 
and  the  complementary  error  function  is  defined  as 


erfc  (x) 


(v2) 


dv  . 


(3) 


Figure  1.  Correlation  detector. 


Denoting  the  inverse  complementary  error  function  by  erfc-1,  we 
define 


6  =  erfc-1  (2Pp)  ,  (4) 

5  =  erfc-1  (2Pd)  •  (5) 

From  equations  (1)  and  (2),  we  can  calculate  the  value  of  E/NQ  necessary 
to  ensure  specified  Pp  and  PD.  The  result  is 

f-  =  (3  -  £)2  •  (6) 

0 

Although  the  ideal  correlator  cannot  be  used  when  s(t)  is  unknown, 
equation  (6)  provides  a  basis  of  comparison  for  more  realistic  intercep¬ 
tion  receivers. 

To  detect  the  presence  of  an  unknown  signal,  we  assume  that  the 
intercepted  signal  has  random  phase  and  frequency  and  an  unknown  con¬ 
stant  amplitude.  The  signal  frequency  is  assumed  to  be  one  of  M  pos¬ 
sible  values;  that  is,  the  band  to  be  searched  is  divided  into  M 
channels  with  center  frequencies  a^,  a>2,  etc.  To  each  discrete  fre¬ 
quency,  ok,  we  assign  a  hypothesis,  Hi>  Thus,  the  multiple  alternative 
hypotheses  over  an  observation  interval  are 
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where  the  9^  are  phase  angles.  We  assume  that  the  phase  angles  are 
uniformly  distributed  and  that  each  frequency  is  equally  likely  to 
occur.  A  comparison  of  the  likelihood  ratios1  yields  the  receiver 
depicted  in  figure  2.  The  decision  rule  is  the  following:  choose  H^, 
i  =  1,  ...»  M,  if  the  largest  envelope  detector  output  is  greater  than 
the  threshold,  and  choose  otherwise.  If  a  signal  is  detected,  this 
receiver  automatically  identifies  the  frequency  as  the  center  frequency 
of  the  filter  with  the  largest  output. 


The  matched  filters  of  figure  2  are  matched  to  intercepted  signals 
that  are  pulsed  sinusoids.  To  accommodate  more  general,  unknown  sig¬ 
nals,  the  matched  filters  could  be  replaced  by  bandpass  filters.  How¬ 
ever,  such  a  replacement  would  give  a  detector  that  is  not  necessarily 
optimum. 


MATCHED 
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Figure  2.  Optimum  detector  for  pulsed  sinusoid  of 
unknown  frequency. 


M.  Whalen,  Detection  of  Signals  in  Noise,  Academic  Press,  Inc.,  New 
York  (1971). 


the  band  |f|  <  W/2.  The  Gaussian  noise  emerging  from  the  bandpass 

filter  can  be  represented  in  terms  of  quadrature  components  as4 

n(t)  =  nc(t)  cos  i*>ct  -  ng(t)  sin  wct  .  (8) 

If  n(t)  is  filtered  white  noise  of  spectral  density  Nq/2,  then  nc(t)  and 
ng(t)  have  flat  power  spectral  densities,  each  equal  to  NQ  over  |f|  < 

W/2. 

As  shown  in  figure  3,  the  input  to  the  comparator  is 

V(t)  =  J^_T  y2(T)  dx  ,  (9) 

where  the  integration  interval  is  equal  to  the  observation  interval,  T. 

• 

The  comparator  output  may  be  sampled  or  continuously  fed  to  a 
processor.  We  next  determine  the  probabilities  of  false  alarm  and 

detection  associated  with  V(t)  at  a  fixed  time.  For  convenience,  we 

normalize  the  test  statistic  to 

V  =  /  y2(t)  dt  .  MO) 

0 

Substituting  equations  (7)  and  (8)  and  assuming  that  fcT  >>  1,  fc  >>  W, 
we  obtain  the  approximation 

T 

V  =  —  /  [s  (t)  +  n  (t)l2  dt 

N  o  L  c  c  J 


+  J  [s  (t)  +  n  (t)l2  dt  .  (ID 

N  o  L  s  s  J 

0 

From  the  ^sampling  theorems  for  deterministic  and  stochastic  processes,5 
respectively,  we  obtain  expressions  that  facilitate  a  statistical  per¬ 
formance  analysis.  After  an  appropriate  choice  of  time  origin,  we  may 
write 

4R.  E.  Ziemer  and  W.  H.  Tranter,  Systems,  Modulation  and  Noise, 
Houghton  Mifflin  Co.,  New  York  (1976). 

5A.  Papoulis,  Signal  Analysis,  McGraw-Hill  Book  Co.,  New  York  (1977). 
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00 


sc(t)  =  l  sc  j  sine  (Wt  -  i)  , 
00 

Sg(t)  =  l  ^  s sine  (Wt  -  i)  , 
00 

nc ( t)  =  l  nc{£)  sine  (Wt  -  i)  , 
00 

ns(t)  =  J  ns(w)  sinc  (Wt  ~  i>  ' 


where  sinc  x  -  sin  ttx/itx.  We  make  the  following  approximations,  based 
upon  the  known  properties  of  the  sinc  function: 


/Q  sinc  (Wt  -  i)  sinc  (Wt  -  j)  dt  =  0  ,  i  *  j  , 


/  sinc2  (Wt  -  i)  dt  -  J°°  sinc2  ( Wt.  -  i)  dt  =  -  ,  0  <  i  <  TW  ,  (17) 

0  W 


/  sinc2  (Wt  -  i)  dt  -  0  ,  i  <  0  or  i  >  TW  .  (18) 

The  error  introduced  by  each  integral  approximation  is  bounded  by 
1/2W.  Assuming  that  TW  >_  1,  the  error  introduced  by  equation  (18)  at 
i  =  0  is  nearly  1/2W.  For  other  values  of  i,  except  possibly  i  =  TW, 
the  errors  caused  by  the  approximations  are  much  less  than  1/2W  and 
decrease  as  TW  increases.  Substituting  equations  (12)  to  (18)  into 
equation  (11),  we  obtain 


V  j,  [*«  (»)  +  *0  (»)]2 

*  *7  ®  *  ”■  (»)  2  ' 


% 


where  Y  is  the  largest  integer  less  than  or  equal  to  TW.  In  view  of  the 
approximations  made,  this  equation  becomes  an  increasingly  accurate 
approximation  of  equation  (11)  as  Y  increases.  It  is  always  assumed 
that  Y  1* 

We  assume  that  the  bandpass  filter  has  a  transfer  function  that 
is  rectangular  about  fc.  Since  n(t)  has  a  power  spectral  density  that 
is  symmetrical  about  f c ,  nc(t)  and  ng(t)  are  independent  Gaussian  pro¬ 
cesses.4  Thus,  nc(i/W)  and  ng(j/W)  are  independent  Gaussian  random 
variables.  The  power  spectral  densities  of  both  nc(t)  and  ng(t)  are 
S(f)  =  Nq  for  |f  |  <  W/2  and  S(f)  =  0  otherwise.  The  associated  auto¬ 

correlation  function  is 


R(T)  =  NqW  sine  WT  .  (20) 

This  expression  indicates  that  nc(i/W)  is  statistically  independent  of 
nc(j/W),  i  /  j,  and  similarly  for  ng(i/W)  and  ng(j/W).  If  n(  t)  is 
assumed  to  be  zero-mean,  so  are  nc( i/W)  and  ng(i/W).  Using  these  facts, 
we  rewrite  equation  (19)  as 

V  =  I  aj  +  l  b2  ,  (21) 

i=1  i=1 


where  the  a^'s  and  the  b^'s  are  statistically  independent  Gaussian 
random  variables  with  unit  variances  and  means 


m1i  =  E  [ai]  = 


(N 0wy/2  c  w 
21  =  Elbi^ =  (^p2  Ss  (*)  * 


(22) 


(23) 


The  first  sum  in  equation  (21)  has  a  noncentral  X2  distribu¬ 
tion1  with  Y  degrees  of  freedom  and  a  noncentral  parameter  XJ  =  ^  m2^. 
Similarly,  the  second  sum  has  a  noncentral  x2  distribution  with  Y 
degrees  of  freedom  and  a  noncentral  parameter  X  =  mfi-  Since  the  two 


lA.  Whalen,  Detection  of  Signals  in  Noise,  Academic  Press,  Inc.,  New 
York  (1971). 

4R.  E.  Ziemer  and  w.  H .  Tranter ,  Systems,  Modulation  and  Noise, 
Houghton  Mifflin  Co.,  New  York  (1976). 
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2 

X  variables  are  independent,  V  has  a  noncentral  X2  distribution  with  2Y 
degrees  of  freedom  and  noncentral  parameter  A  =  A^  +  A2.  Thus, 


A  = 


■  ;r  K  [•=“’  *  •J<t>]  at 

2  T 

~  —  /  s2(  t)  dt  . 

No  ° 


(24) 


In  terms  of  the  signal  energy,  E,  we  have  the  approximation 


0 


(25) 


By  straightforward  calculations  using  the  statistics  of  Gaussian  var¬ 
iates,  the  mean  and  the  variance  of  V  are  determined  to  be 


E  [V]  =  A  +  2Y  ,  (26) 

VAR(V)  *  4A  +  4Y  .  (27) 

By  using  the  known  probability  density  functions  for  a  noncentral  X2 
random  variable,  the  false  alarm  and  detection  probabilities  can  be 
expressed  as  integrals.  In  the  absence  of  a  signal,  the  X2  probability 
density  function  for  V  is 

p.(v)  =  -y-J —  vY_1  e“v/2  ,  v  i  0  (28) 

0  2Tr(Y) 

P0(v)  =  0  ,  v  <  0 


where  T(x)  is  the  gamma  function.  The  false  alarm  probability  is 


P q  ( v)  dv  . 


(29) 


12 


.1 


1 1 


* ' 


If  the  signal 

v) 

p  (v) 
1 


is  present. 


the  X2 


1Q(Y-i)/2 


exp 


=  0  , 


probability  density  function  for  V  is 

(-  H1)  v, </Wl  •  v  >  o . 

v  <  0  ,  oo) 


where  In(x)  is  the  modified  Bessel  function  of  the  first  kind  and  order 
n.  The  probability  of  detection  is 


P D  =  /  P  j  ( v )  dv 


(31) 


Numerous  schemes  for  evaluating  Pp  and  PD  have  been  proposed  in  the 
literature . * ' 3 

We  are  particularly  interested  in  the  case  in  which  TW  is  large 
since  this  case  includes  the  interception  of  spread  spectrum  communica¬ 
tions.  When  TW  is  large,  Y  -  TW,  and  the  central  limit  theorem  indi¬ 
cates  that  V  is  approximated  by  a  Gaussian  variate.  Using  equations 
(26)  and  (27)  with  A  =  0  and  equation  (29)  with  a  Gaussian  density  for 
Pq(v) ,  we  obtain 


F  (8htwJ1//2 


exp 


(v  ~  -2.TW) 

omr.7  OV 


8TW 


erfc 


VT  -  2TW 
(8TW)  lj/2 


(32) 


Similarly,  we  obtain 


PD  =  2  erfc 


VT  -  2TW  -  A 


(8TW  +  8A) 


172 


(33) 


1A.  Whalen,  Defection  of  Signals  in  Noise,  Academic  Press,  Inc.,  New 

York  (1971). 

3tf.  Urkowitz,  Energy  Detection  of  Unknown  Deterministic  Signals, 
Proc .  IEEE,  55  (April  1967),  523. 
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Combining  equations  (4),  (5),  (32),  and  (33)  gives 


( 8TW  +  8X)1/2  £  =  (8TWJ1/2  8  -  A  . 


(34) 


We  can  solve  this  equation  to  determine  the  value  of  E/Nq  necessary  to 
achieve  specified  values  of  Pp  and  PD»  Solving  for  A  and  using  equation 
(25),  we  obtain 


E_ 
N  r 


r  l1'2 

=  2C2  +  8(2TW)1/>2  -  s|_2TW  +  4€2  +  23(8TW)1/2J 


(35) 


If  we  assume  that  TW  >>  82  and  TW  >>  C2,  then  further  simplification  is 
possible.  The  result  is 

f-  =  (2TW)1/2  (3  -  5)  ,  TW  »  max  (32,  C2}  .  (36) 

N0 

Comparing  equations  (36)  and  (6),  we  see  that  the  disparity  in  perform¬ 
ance  between  the  radiometer  and  the  matched  filter  increases  with  TW. 
Equation  (36)  indicates  that  detection  difficulties  increase  as  the 
intercepted  signal  spectum  is  spread. 

Denoting  the  intercepted  signal  power  by  Rg  and  the  signal  duration 
by  Tj,  the  intercepted  power  necessary  to  achieve  specified  values  of  Pp 
and  PD  is 

Rs  =  N0  (3  -  ?)  ,  T  L  <  T  ,  TW  »  max  ( 82,  Z2)  , 

Rs  =  No(^t)1/2  "  5)  r  T1  >  T  ,  TW  »  max  ( 32  ,  ?2)  .  (37) 


As  long  as  T ^  T,  this  equation  indicates  that  increasing  the  observa¬ 
tion  interval  decreases  the  required  power.  However,  if  T^  <  T,  an 
increase  in  the  observation  interval  increases  the  required  power. 

If  the  outputs  of  Nr  independent  radiometers  are  averaged,  a 
straightforward  calculation  shows  that  the  required  Rg  can  be  reduced  by 
a  factor  of  Nr~1//2. 


2.2  Channelized  Radiometer 


A  channelized  radiometer  forms  when  M  radiometers  are  inserted 
in  the  branches  of  figure  2,  as  depicted  in  figure  4(a).  Each  block 
labeled  radiometer  contains  a  bandpass  filter  of  bandwidth  W/M,  a  squar¬ 
ing  device,  and  an  integrator,  but  no  comparator.  Let  Tg  denote  the 
sampling  interval,  which  is  the  observation  interval  of  the  constituent 
radiometers.  To  avoid  processing  extraneous  noise,  the  arrival  time  of 
the  signal  to  be  intercepted  may  be  estimated  by  additional  hardware. 
The  sampling  interval  may  equal  or  be  somewhat  less  than  the  minimum 
expected  signal  duration  in  a  channel .  To  increase  effectiveness 
against  frequency  hopping  or  multiple  frequency-shift  keying  (MFSK),  the 
processor  examines  N  consecutive  comparator  outputs  and  determines  that 
a  signal  is  present  if  r  of  these  outputs  correspond  to  comparator 
inputs  that  exceed  the  threshold.  For  example,  if  N  is  odd,  a  majority 
decision  rule  requires  r  =  (N  +  1)/2.  The  effective  observation  inter¬ 
val  of  the  channelized  radiometer,  given  by  T  =  NTg ,  should  usually  be 
less  than  the  minimum  expected  message  duration.  If  it  is  known  that 
the  intercepted  signal  is  narrowband,  we  can  set  N  =  1.  If  the  presence 
of  more  than  one  signal  is  to  be  verified,  it  is  desirable  to  employ  an 
array  of  radiometers  of  the  form  of  figure  3  with  the  comparator  outputs 
feeding  into  a  processor  that  analyzes  the  activity  of  individual  chan¬ 
nels,  as  shown  in  figure  4(b).  In  this  configuration,  N  =  1  and  each 
bandpass  filter  has  a  bandwidth  of  W/M. 


i.i 


lb) 


Figure  4.  Channelized  radiometers  (a)  for 
detection  of  presence  of  hostile  communi¬ 
cations  and  (b)  for  simultaneous  detection 
of  multiple  signals. 
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To  simplify  the  analysis  of  the  interception  of  a  single  sig¬ 
nal,  we  assume  that  the  N  sets  of  radiometer  outputs  are  statistically 
independent.  If  Pp  is  the  probability  that  a  particular  radiometer 
output  exceeds  the  threshold  when  no  signal  is  present,  then  the  proba¬ 
bility  that  none  of  the  radiometer  outputs  exceeds  the  threshold  is  (1  - 
Pp)M,  assuming  that  the  channel  noises  are  statistically  indepen¬ 
dent.  The  probability  that  exactly  i  out  of  N  comparator  inputs  exceeds 
the  threshold  is 


p  u.  »  -  (?)  [i  -  (i  -  p')  -  p’) 


M(N-i) 


i  <  N 


P  (i,  N)  =  0  , 


i  >  N  . 


It  follows  that  the  probability  of  false  alarm  associated  with  the 
observation  interval  is 


Pp  -  l  P(i,  N)  . 
i=r 


If  the  intercepted  signal  duration,  T^ ,  is  less  than  the  observation 
interval,  T,  we  assume  for  simplicity  that  N1  =  Tj/Tg,  the  number  of 
sampling  intervals  during  which  the  signal  is  present,  is  an  integer. 
Further,  we  assume  that  a  single  radiometer  contains  the  intercepted 
signal  during  each  sampling  interval.  he t  P£  denote  the  probability 
that  the  threshold  is  exceeded  at  the  end  of  a  sampling  interval  when  a 
signal  is  present.  Let  P£>  denote  the  probability  that  a  particular 
radiometer  output  exceeds  the  threshold  when  a  signal  is  present  in  that 
radiometer.  From  these  definitions,  it  follows  that 


ps- 1  -(’  -  py  c  -*■?)' 


The  probability  of  detection  associated  with  the  observation  interval  is 
determined  by  reasoning  similar  to  that  which  led  to  equations  (38)  and 
(39).  The  result  is 


?D  =  ^  0  -  Pd)  1  3  P(i  -  j'  N  -  \)  • 

i=r  j=0  \ 


To  compare  the  channelized  radiometer  with  a  single  wideband 
radiometer,  we  assume  that  N  =  N1  and  that  the  energy  of  the  intercepted 
signal  is  equally  divided  among  the  N  sampling  intervals.  Since  the 
total  receiver  bandwidth  is  W,  the  bandwidth  of  each  constituent  radi¬ 
ometer  is  Wg  =  W/M.  Thus,  for  large  values  of  TgWg,  P£  and  P^  are  given 
by  equations  (32)  and  (33)  with  Wg  =  W/M  substituted  for  W,  T  =  T/N 
substituted  for  T,  and  Ag  =  A/n  substituted  for  A.  we  have 


.  MNV_  -  2TW 

P£  -  f  erfc  - ,  (42) 

^  (8TWMN)i/<d 


We  de  f ine 

31  =  erfc-1  (2pf)  '  (44) 


C1  =  erfc-1  (2pd)  •  (45 ) 


If  Pp  and  PD  are  specified,  we  solve  equations  (38)  to  (41)  for 
Pp  and  P^.  Using  equations  (42)  to  (45),  we  perform  a  calculation 
analogous  to  that  used  in  deriving  equation  (37).  The  result  is  that 
the  required  Rg  for  detection  with  specified  values  of  Pp  and  PD  is 

Rs  a  No  (^mt)  (ei  -  ^l)  '  T1  ±  T,  TW  »  MN  max  ^  .  (46) 

If  a  frequency-hopping  signal  is  to  be  intercepted,  the  parameter  N  is 
proportional  to  the  hopping  rate.  Thus,  equation  (46)  indicates  that 
the  required  power  is  proportional  to  the  square  root  of  the  hopping 
rate.  If  M  »  N,  the  channelized  radiometer  requires  less  power  than  a 

wideband  radiometer  with  the  same  values  of  T  and  W.  If  we  set  N  =  1 

and  attempt  to  intercept  a  frequency-hopping  signal  by  processing  each 
hop,  then  T  must  be  decreased  as  the  hopping  rate  increases. 
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Suppose  that  the  energy  is  concentrated  in  a  narrow  bandwidth 
during  the  observation  interval  and  that  the  bandwidths  of  the  radi¬ 
ometers  are  sufficiently  wide  so  that  the  energy  enters  a  single  radiom¬ 
eter.  Then  we  may  take  r=N=N^=1.  Equations  (39)  and  (41)  reduce 
to 


The  required  value  of  Rg  for  detection  is  determined  by  the  usual  method 
to  be 


T,  >  T,  TW  »  M  max  (3 


Thus,  the  required  power  falls  approximately  as  the  square  root  of  the 
number  of  channels. 


To  determine  P£,  and  when  equations  (42)  and  (43)  do  not 
apply,  we  assume  that  the  intercepted  signal  energy  is  equally  divided 
among  the  Nj  sampling  intervals.  We  use  equations  (28)  to  (31)  with  Xg 
=  X/Nj  substituted  for  X  and  h,  the  largest  integer  less  than  or  equal 
ho  T  W  ,  substituted  for  Y.  The  results  are 


By  numerical  methods,  this  pair  of  equations  can  be  solved  simultan¬ 
eously  to  eliminate  VT  and  express  A  =  NjAa  as  a  function  of 
and  Pjj  .  If  we  solve  equations  (38)  to  (41)  for  P^  and  P^  in  terms  of 
PF  and  PD,  then  we  can  obtain  an  equation  for  A  in^terms  of  PF  and  PD> 
From  this  equation,  we  finally  obtain  the  required  value  of  Rg  necessary 
to  achieve  specified  values  of  Pp  and  PD< 

The  channelized  radiometer  has  been  shown  to  be  relatively 
effective  against  conventional  and  frequency-hopping  communications.  It 
is  also  useful  against  pseudonoise  spread  spectrum  communications  if 
preliminary  processing  is  used  to  produce  a  signal  with  a  narrow  band¬ 
width  (sect.  3). 

2.3  Cross  Correlator 


The  ideal  correlator  of  figure  1  can  be  approximated  if  the 
signal  is  intercepted  at  two  spatially  separated  antennas.  The  cross¬ 
correlation  function  of  the  two  antenna  outputs  is  estimated  for  various 
relative  arrival  times,  and  the  peak  value  of  this  function  is  applied 
to  a  comparator.  Figure  5  is  a  block  diagram  of  a  realization  employing 
the  discrete  Fourier  transform  (DFT).  One  way  to  implement  the  DFT  is 
to  use  a  digital  filter  and  the  fast  Fourier  transform  algorithm.  An 
alternative  implementation  is  to  use  the  chirp  Z- transform  algorithm  and 
charge-coupled  devices.  An  analog  realization  of  the  cross-correlation 
function  that  is  similar  to  the  configuration  of  figure  5  can  be  accom¬ 
plished  with  dispersive  filters  providing  Fourier  transforms  (sect. 
3.6).  Elegant  realizations  are  possible  with  acousto-optical  devices.6 

Figure  5a  depicts  the  initial  processing  of  each  antenna  out¬ 
put.  After  passage  through  the  bandpass  filter  of  bandwidth  W,  the 
intercepted  waveform,  r(t)  =  s(t)  +  n(t),  can  be  represented  as 


r(t)  =  rc(t)  cos  wct  -  rg(t)  sin  wct 


(54) 


where  the  quadrature  components,  rc(t)  and  rg(t),  are  confined  to  the 
band  |f|  <_  W/2.  If  the  two  lowpass  filters  have  bandwidths  W/2  and  f 
=  w c/2 tt  >  W/2,  then  rc(t)  and  rg(t)  are  extracted  by  the  operations 
shown  in  figure  5(a).  Analog-to-digital  converters  produce  the  discrete 
sequences  rc( i/W)  and  rg(i/W). 

6F.  A.  Sprague,  A  Review  of  Acousto-optic  Signal  Correlators ,  Optical 
Engineering ,  16_  (September  1977),  467. 
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-2  sin  <uct 

(a) 


<b) 


Figure  5.  Cross  correlator:  (a)  initial  processing  of  each  antenna 
output  and  (b)  final  processing. 


We  regard  one  of  the  antennas  as  a  reference  and  denote  its 


output  by  the 
have 

subscript  1. 

In  terms 

of  the 

signal 

and  the  noise,  we 

r  j  ( t) 

=  s(t)  + 

nj  (t) 

9 

(55) 

r1c(t) 

=  sc(t)  + 

n1c(t) 

9 

(56) 

=  ss(t)  + 

n1s(t> 

• 

(57) 
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r2(t)  =  s(t  -  Tj 

1+  n^(t)  , 

(58) 

2c<t)  =  sc(t  -  Tr] 

|  +  n2c(t)  ' 

(59) 

2s(t)  =  sg(t  -  Tr) 

|  +  n2g( t)  , 

(60) 

where  Tr  is  the  arrival  time  of  the  intercepted  signal  at  this  antenna 
output  relative  to  the  arrival  time  at  the  reference  antenna  output.  By 
inserting  a  sufficiently  long  delay  before  this  antenna  output,  we 
ensure  that  T_  >  0 .  It  is  convenient  to  use  the  notation 


i  =  0,  1,  .  .  Y  -  1 


for  two  of  the  discrete  sequences  observed  during  the  interval  T.  We 
form  sequences  with  K  =  y  +  Nc  -  1  points  by  augmenting  the  a^  and  b^ 
1  zeros,  where  N  <  y .  As  indicated  in  figure  5(b)  ,  the 


conjugate  DFT  of  a^  is  calculated,  giving 


K-1 

v  -in 

A?  =  >  a 

*  n  n  K 

n=0 


where  (2  =  exp  (-j2ir/K)  and  j  =  /-I  .  Similarly,  the  DFT  of  b^  is 


\  *  \  bnaKn  ■  1  '  ’ .  K  -  ’  ' 

n=0 
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(63) 


The  inverse  DPT  of  the  product  A*B^  for  0  <  Nc  -  1  is 

.  K-1 

C.  =  -  l  A*B  Q  ~ln 
1  K  n  n  K 


it  l„  X  aA  *, 

m=0  k=0  n=0 


Y_1  y_i  Y-1 

=  i  y  a  bT,^  t  • =  y  a  .  b  , 

„  'v  m  k  m,k-i  ,  ^ .  k-i  k 

k=0  m=0  k=i 


where  6^^  is  the  Kronecker  delta.  From  the  original  definitions  and 
equations  similar  to  equations  (12)  to  (18),  we  obtain  for  i  <<  TW, 


ci  =  J  rlc(k  H  "■■)r2c(HJ) 

T 

“"{/»r,^t"”)r2<=,t>  at  '  1651 

This  sequence  is  the  output  of  one  of  the  inverse  DFT  operations  show.i 
in  figure  5(b) .  An  analogous  expression  can  be  written  for  the  output 
of  the  other  inverse  DFT  operation.  The  addition  of  the  two  sequences 
produces  a  sequence  proportional  to 

Di  *  {,  M*  -  iKc(t>  +  r1e(‘  -  iK'11  dt  '  <66’ 

i/W  l_  ' 


If  f  >>  W,  then  expansions  similar  to  equation  (54)  lead  to  the  approx¬ 
imation 


D  =  /  2r  ft  -  -7)  r  (t)  dt  ,  i  =  0,  1,  .  .  .,  N  -  1 

1  4  1  V  2  c 
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To  interpret  the  next  operation  in  figure,  5(b),  we  initially 
assume  that  no  noise  is  present.  In  this  case,  equations  (55)  and  (58) 
yield 


T 

D  =/  2s  t  -  i  s  t  -  T  dt  ,  i  =  0,  1,  .  .  .,  N  -1  .  (68) 

i/W  w  r  C 

Thus,  the  Di  provide  sampled  values  of  an  approximation  of  the  auto¬ 
correlation  function  of  s(t).  Let  iQ  denote  the  index  that  corresponds 
to  the  largest  .  Assuming  that  the  approximation  is  adequate  and  that 
the  autocorrelation  function  has  a  sharp  peak,  i  is  the  index  closest 
to  the  value  TrW.  When  noise  is  present,  this  statement  may  not  be 
true;  however,  to  proceed  with  the  analysis,  we  assume  that  it  is.  Note 
that  Tr  can  be  estimated  as  in/W.  This  estimate  can  be  used  for 
direction  finding  (sect.  4.2).  u 

Assuming  that  the  largest  Di  has  index  i  =  T  W  <  N  -  1  and 
normalizing,  the  input  to  the  comparator  in  figure  5(b)~isCthe  test 
statistic 


V  =  £  /  r  (t  -  T  )  r  (t)  dt 

N„  m  1V  ^2 


0  T 


r 
T-T 


=  —  /  r  r  (t)r  (t  +  T 
*0  0  1  r) 


dt  . 


(69) 


"N 


Substituting  equations  (55)  and  (58)  and  defining  n2(t)  =  n'(t  +  Tr) 
and  Ta  =  T  -  Tr,  we  get  2 


V  =  j  a  [s(t)  +  n  (t)l  [s(t)  +  n  (t)l  dt 
Nq  o  L  !  J  L  2  J 


(70) 


In  the  usual  manner,  we  obtain  the  series  expansion 


0  i=1 


(w)  +  nic  sc(^)  +  n2c(wj 


—  I 

NoW  i=1 


(w)  +  nisfe)  Ss6)  +  n2s(i) 


+ 


(71) 


where  Ya  is  the  largest  integer  less  than  or  equal  to  TaW.  Assuming 
that  n^t)  and  n2(t)  are  statistically  independent,  zero-mean,  Gaussian 
processes,  a  straightforward  calculation  yields 


E[V]  =  A  , 

a 


VAR(V)  =  2A  +  2Y  , 

cl  3 


where  A  =  =  2E  /Nn  and  Ea  is  the  energy  in  interval  T  . 

d  d  U  u  a 

For  large  values  of  TflW,  the  test  statistic  is  approximately 
normally  distributed.  It  follows  that 


pD  =  j  erfc 


Vm  -  A 
T  a 


(4TaW  + 


4*a) 


171 


(74) 


According  to  equation  (67),  the  cross-correlation  function  is 

However,  the  maxim  on  possible  value  of 


computed  for  Nc  sample  values. 


Tr  way  be  such  that  only  Nc  <<  y  sample  values  need  to  be  computed  to 
obtain  V.  When  no  signal  is  present,  a  false  alarm  occurs  if  any  of  the 
Nc  estimated  cross-correlation  values  exceeds  the  threshold.  If  Nr  is 
sufficiently  small,  it  is  reasonable  to  assume  that  each  estimated  value 
has  approximately  the  same  probability,  denoted  by  P^,  of  exceeding  the 
threshold.  This  assumption  implies  the  approximation. 


Equations  (44)  and  (75)  yield 


=  erfc-1  [2  -  2(l  -  Pp)1/Nc] 


(78) 


We  obtain  in  the  usual  manner  the  required  Rg  to  detect  a  signal 
with  specified  values  of  Pp  and  PD.  The  result  is 

(t  /  \  / 

N0  (B1  “  7  '  T1  <  Ta'  TaW  >:>  max  (61'  ' 

(79) 

=  N°(^)1/2  (3l  "  *)'  T1  2  Ta  '  TaW  »  (Bl'  «*)  ' 

where  Tj  is  the  signal  duration. 

Comparison  with  equation  (37)  indicates  that  the  cross  corre¬ 
lator  can  give  a  theoretical  improvement  of  approximately  1.5  dB  over  a 
single  wideband  radiometer .  Taking  into  account  the  approximations  made 
to  derive  equation  (79),  it  is  possible  that  in  practice  the  cross 
correlator  provides  no  improvement  at  all.  A  comparison  of  figures  3 
and  5  indicates  that  the  implementation  of  the  cross  correlator  entails 
considerably  more  hardware  than  the  implementation  of  a  wideband  radi¬ 
ometer.  However,  as  discussed  in  subsequent  sections,  the  cross  corre¬ 
lator  requires  little  additional  hardware  to  provide  frequency  esti¬ 
mation  and  direction  finding. 

The  channelized  cross  correlator  is  an  array  of  M  cross  corre¬ 
lators,  each  of  which  has  a  bandwidth  of  W/M.  The  outputs  of  the  array 
are  applied  to  a  processor.  Analogously  to  the  channelized  radiometer, 
the  channelized  cross  correlator  may  be  preferable  to  a  single  wideband 
cross  correlator  when  the  hostile  communications  are  narrowband  or  when 
tv«  or  more  simultaneous  signals  are  to  be  intercepted. 

Equations  (37)  and  (79)  indicate  that  increasing  the  bandwidth 
of  a  frequency-hopping  system  degrades  the  performance  against  a  single 
signal  of  both  the  wideband  cross  correlator  and  the  wideband  radiom¬ 
eter.  However,  neither  of  these  receivers  is  sensitive  to  the  hopping 


R 
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rate.  Increasing  the  hopping  rate  makes  the  practical  design  of  a 
channelized  receiver  more  difficult  and  degrades  its  performance.  If 
the  rate  is  sufficiently  high,  the  channelized  receiver  may  have  to  be 
abandoned  in  favor  of  a  wideband  receiver. 


3.  FREQUENCY  ESTIMATION 

The  immediate  purpose  of  a  frequency  estimation  system  is  to  deter¬ 
mine  the  center  frequency  and  possibly  the  spectral  shape  of  an  inter¬ 
cepted  signal.  If  a  frequency-hopping  signal  is  intercepted,  the 
purpose  is  to  determine  each  hopping  frequency  or  at  least  the  frequency 
range  over  which  the  hopping  occurs. 

Although  not  desirable  for  some  purposes,  such  as  message  analysis, 
Preliminary  processing  of  pseudonoise  spread-spectrum  communications  is 
desirable  before  estimation  of  the  center  frequency  is  attempted.  An 
intercepted  binary  pseudonoise  signal  has  the  form 

s(t)  =  Am(t)p(t)  cos  w  t  ,  (80) 

0 

where  A  is  the  amplitude,  is  the  center  frequency,  m(t)  is  the 

binary  message  sequence,  and  p(t)  is  a  binary  pseudorandom  sequence. 
Both  m(t)  and  pit )  take  the  values  +1  or  -1.  Suppose  s(t)  enters  a 
wideband  receiver  and  is  squared.  Since  m2  =  p2  =  1 ,  the  output  of  the 
squaring  device  is  proportional  to 

2  2 

s2(t)  =  A2  cos2  w  t  =  cos  2  t  .  (81) 

0  2  2  0 

The  double- frequency  term  is  now  a  pure  pulsed  sinusoid.  Its  frequency 
and,  thus,  the  center  frequency  of  s(t)  can  be  estimated  by  the  systems 
described  in  this  section.  The  same  preliminary  processing  is  useful 
against  phase-shift  keyed  communications . 

3. 1  Channelized  Receiver 

Estimation  theory  leads  to  the  receiver  of  figure  2  for  fre¬ 
quency  estimation,  assuming  that  the  arrival  time  and  the  signal  wave¬ 
form,  except  for  a  uniformly  distributed  phase  angle,  are  known.1  After 
the  largest  output  is  selected,  the  unknown  frequency  is  estimated  as 
the  center  frequency  of  the  filter  producing  the  largest  output.  A 
practical  approximation  to  this  receiver  is  the  channelized  radiometer 
of  figure  4. 

ZA.  Whalen,  Detection  of  Signals  in  Noise,  Academic  Press,  Inc.,  New 
York  (1971). 
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Suppose  we  desire  a  frequency  resolution  of  A,  where  A  is  not 
less  than  the  Cramer-Rao  bound.1  If  the  entire  range  of  reconnoitered 
frequencies  is  W,  then  each  filter  must  have  bandwidth  2A  and  M  =  W/2A 
is  the  number  of  required  filters  to  attain  the  desired  accuracy.  If 
the  intercepted  signal  has  duration  T,  then  each  filter  must  have  band¬ 
width  2/T  for  most  of  the  signal  energy  to  pass  through  it.  Thus,  we 
have 


(82) 


If  M  >_  6,  the  number  of  required  filters  can  be  reduced  by 
arranging  the  filters  in  successive  stages,  as  shown  in  figure  6  for  the 
case  in  which  each  stage  has  the  same  number  of  filters.  After  an 
intercepted  signal  passes  through  the  first  filter  bank,  its  frequency 
is  theoretically  known  within  an  accuracy  of  A^  =  w/2Mg,  where  Mg  is  the 
number  of  filters  in  each  stage  and  2A^  is  the  bandwidth  of  the  first 
stage  filters.  A  bank  of  mixers  ensures  that  the  filter  outputs  are 
shifted  in  frequency  so  that  the  input  to  the  second  stage  has  a  fre¬ 
quency  between  fn  -  fc1  -  A  and  f^  -  fc1  +  Ajf  where  f  is  the 
center  frequency  of  the  top  filter  in  the  first  bank  and  f^  is  the 
frequency  of  a  local  oscillator.  After  the  input  passes  through  the 
second  filter  bank,  the  frequency  is  known  within  an  accuracy  A^  = 
A^/Mg,  and  so  on.  If  Ng  stages  of  Mg  filters  each  are  employed,  then  an 
accuracy  of  A  is  attained  if 
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M 


W_ 

2A 


(83) 


The  total  number  of  filters  required  is  MgNg.  Disadvantages  of  the 
channelized  receiver  of  figure  6  relative  to  that  of  figure  4  are  the 
increased  processing  time  required  for  frequency  estimation,  the  reduced 
amount  of  noise  and  interference  filtering,  and  the  ambiguities  that 
arise  when  more  than  one  signal  is  intercepted. 


It  is  not  necessary  that  each  stage  have  the  same  number  of 
filters.  However,  if  Ng,  W,  and  A  are  fixed,  it  can  easily  be  shown, 
using  Lagrange  multipliers,  that  the  total  number  of  filters  is  mini¬ 
mized  if  each  stage  has  approximately  the  same  number  of  filters  (exact¬ 
ly  the  same  number  if  an  integer  Mg  exists  that  satisfies  eq  [83] )  .  If 
each  stage  has  the  same  number  of  filters  and  W  and  A  are  fixed,  it  can 
be  shown  that  the  total  number  of  filters  is  minimized  if  each  stage  has 
three  filters.  (Lagrange  multipliers  yield  Mg  =  e,  but  Mg  must  be  an 
integer. ) 


lA.  Whalen,  Detection  of  Signals  in  Noise,  Academic  Press,  Inc.,  New 
York  (1971). 
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Figure  6.  Channelized  receiver  with  filters  arranged  in 
successive  stages. 

Even  when  minimized,  the  number  of  filters  and  detectors 
required  in  a  channelized  receiver  may  make  this  method  of  frequency 
estimation  expensive.  In  this  case,  a  limited  number  of  filters  can  be 
used  to  reduce  the  total  bandwidth  examined  by  other  frequency  estima¬ 
tion  devices. 

3.2  Spectrum  Analysis  with  Discrete  Fourier  Transform 


The  outputs  of  the  two  conjugate  DFT  blocks  in  figure  5(b) 
provide  a  scaled  phase-shifted  estimate  of  the  intercepted  spectrum  over 
the  receiver  bandwidth.  To  see  how  the  spectrum  is  estimated,  we  first 
define  the  truncated  waveforms: 


s' ( t)  =  s(t)q(t)  , 


sMt)  =  sc(t)q(t)  , 


s' (t)  =  s  (t)q(t) 


(84) 


where 


q(t)  =  1  ,  0  <  t  _<  T  , 

q(t)  =  0,  otherwise  . 


(85) 


We  denote  the  Fourier  transforms  of  s'(t),  s'(t),  and  s'(t)  by  S(w) , 

S  (<*>)  ,  and  S  (to),  respectively.  From  equation  (7), 
c  s 

S(<Jd)=-^S^(D-U)^+4s  (to  +  0)  ^ 

2  c  V  c  /  2  c  V  cj 

—  ~T  S  (to  -  U)  \  +  ^tS  (us  +  to  \  .  (86) 

2)  s\  c)  2]  si  c) 

We  are  Interested  in  determining  S(w)  for  w  >  0.  If  u)  >  ttW,  then 

—  c  — 

S(“)  =  ls  to  >  0  ,  to  >  trw  .  (87) 

2  c  V  c)  2]  s  \  c)  —  c  — 


Thus,  we  can  estimate  S(w)  if  we  first  estimate  Sc(m) 


and  Sg(w) . 


We  give  the  details  of  the  estimation  of  Sc(u))j  the  estimation 
of  Ss ( oj )  is  similar.  For  simplicity,  we  assume  K  =  TW.  The  sample 
values  of  S  (w)  are  related  to  those  of  s'(t)  through  the  fundamental 

,  .  c  C  C 

relation3 


where  0  =  exp  (~j2F/N)  and 


s  (my)  =  l  s’ (my  +  kNy) 
c  ,  _  c 

k=-°° 


S  (nx)  =  S  (nx  +  iNx) 

c  «  c 


(89) 


(90) 
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We  apply  these  equations  with  N  =  Y  =  TW,  x  =  2VT,  and  y  =  1/W.  Since 
s^(t)  =  0  unless  0  <  t  £  T,  equation  (89)  implies  that 


(91  ) 
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If  no  noise  is  present,^  then  Sc(to)  =  Sc(oj)  at  the  sample  values.  The 
resolution  of  Sc(io)  and  Ss(u)),  the  corresponding  estimator  for  S  (u>)  ,  is 
primarily  determined  by  the  duration  of  the  sample  pulse?  The 
resolution  is 


(96) 


To  obtain  a  more  accurate  expression  for  A,  the  effects  of  the  random 
noise  must  be  evaluated.  However,  equation  (96)  is  adequate  for  roughly 
comparing  the  frequency  estimation  potential  of  the  cross  correlator 
with  that  of  competitive  systems . 

If  the  frequencies  of  frequency  hopping  or  MFSK  communications 
are  to  be  successfully  estimated,  the  observation  time,  T,  must  be  less 
than  the  period  between  frequency  changes. 


3.3  Acousto-optical  Receiver 


Spectrum  analysis  using  acousto-optical  diffraction  has  the 
potential  capability  for  real-time,  wideband  frequency  estimation  of 
many  simultaneous  signals.  The  principal  components  of  an  acousto- 
optical  spectrum  analyzer  are  shown  in  figure  7.  The  diffraction 
geometry  associated  with  the  Bragg  cell  is  illustrated  in  figure  8.  The 
Bragg  cell  converts  an  electronic  input  at  frequency  f  into  a  traveling 
acoustic  wave  with  velocity  vfl  and  wavelength  Aa  =  va/f0.  The  laser 
light  has  wavelength  AQ  in  free  space  and  A  /n  inside  the  cell,  which 
has  an  index  of  refraction  n.  According  to  Bragg's  law,  the  sound  wave 
interacts  with  the  light  beam  to  produce  a  principal  diffracted  beam, 
which  is  offset  from  the  incident  beam  by  an  angle 


0' 


2  sin-1 


(97) 


inside  the  cell  and  an  angle 


0=2  sin-1 


(98) 


outside  the  cell.  These  equations  are  valid  provided  that  the  acoustic 
wave  has  a  single  wavelength  across  the  cell.  For  small  values  of  the 
argument,  equation  (98)  becomes 


A  f 
o  0 

9  -  -  ,  A  f  «  2v  .  (99) 

va  o  o 
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The  lens  produces  a  Fourier  transform  on  its  focal  plane  at  the 
photodetector  array.  The  center  of  the  diffracted  beam  converges  to  a 
position  a  distance 


from  the  center  of  the  corresponding  undiffracted  beam 
focal  length  of  the  transform  lens.  Thus,  the  freq 
estimated  by  measuring  the  relative  intensities  at  t 
array  elements. ' 


:ht,  Spectrum  Analysis  Using  Acousto-optic  Filters,  Optical 
16  (September  1977),  461, 


V: 


The  diffracted  beam  has  an  angular  width  on  the  order  of 
A  /D,  where  n  is  the  effective  aperture  of  the  Bragg  cell. 
Consequently,  the  diffracted  beam  spreads  over  length  of  FA  /D  in  the 
focal  plane.  The  frequency  resolution  is  defined  to  be  the  difference 
in  frequency  between  two  intercepted  signals  such  that  the  corresponding 
positions  in  the  focal  plane  differ  by  the  spread  of  the  diffracted  beam 
in  the  focal  plane.  From  this  definition  and  equation  (100),  the 
resolution  is 


A 


(101) 


where  Tc  is  the  time  that  it  takes  an  acoustic  wave  to  cross  the  cell 
aperture.  This  equation  applies  when  an  acoustic  wave  of  fixed 
wavelength  occupies  the  aperture.  A  necessary  condition  for  its 
validity  most  of  the  time  is  that  Tc  be  less  than  the  period  between 
frequency  changes  of  a  frequency-hopping  signal  or  an  MFSK  signal . 
Thus,  the  resolution  is  no  better  than  the  inverse  of  the  hopping 
period , 


In  frequency  estimation,  only  one  spatial  dimension  of  the 
Bragg  cell  is  used.  It  is  possible  to  design  a  two-dimensional  Bragg 
cell  array  to  estimate  frequency  and  direction  of  arrival  simultan¬ 
eously,  as  shown  in  the  schematic  diagram  of  figure  9.  The  cell  inputs 
are  the  outputs  of  spatially  separated  antennas.  The  intensity  distri¬ 
bution  across  the  photodetector  array  has  one  or  more  maxima  that  are 
vertically  deflected  proportionally  to  the  intercepted  signal  fre¬ 
quencies  and  are  horizontally  deflected  proportionally  to  the  signal's 
angle  of  arrival.  If  the  largest  of  the  outputs  of  the  photodetector 
elements  is  compared  with  a  threshold,  the  presence  of  an  intercepted 
signal  can  be  determined.  Thus,  in  principle,  an  acousto-optical  system 
can  detect  hostile  communications,  estimate  their  frequencies,  and  find 
their  directions. 

3.4  Instantaneous  Frequency  Measurement 

The  instantaneous  frequency  measurement  (IFM)  receiver, 
illustrated  in  figure  10,  is  often  used  to  estimate  radar  frequencies. 
It  is  possible  to  use  it  as  a  supplementary  frequency  estimator  for 
communications,  but  usually  not  by  itself.  Its  operation  is  based  on 
the  relationship  among  carrier  frequency,  path  length,  and  phase  shift 
of  a  signal.  Suppose  that,  after  passage  through  the  bandpass  filter  of 
bandwidth  W,  an  intercepted  signal  has  the  form 

s(t)  =  A(t)  cos  +  <t>(t)J  ,  (102) 

A .  Coppock,  R.  F .  Croce,  w.  L.  Regier,  Bragg  Cell  RF  signal 
Processing ,  Microwave  J.  (September  1978),  62. 
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where  A(  t)  is  the  amplitude  modulation  and  4>(t)  is  the  angle  modulation 
function.  As  shown  in  figure  10,  this  signal  is  delayed  by  time  6  in 
one  branch  relative  to  the  other  branch.  If  6  is  sufficiently  small, 
then  A(t  -  6)  ~  A(t)  and  4>(t  -  6)  =  4>(t)  for  most  of  the  time.  It 
follows  that 


s j ( t)  ~  A(t)  cos 


|wQt  +  <j>(t)  -  wQ6j 


A(t  -  6)  =  A(t) 


<J>(t  -  6)  =  <j>(t)  . 


(103) 
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Figure  9.  Array  of  Bragg  cells  for  simultaneous  frequency 
estimation  and  direction  finding. 


Figure  10.  Instantaneous  frequency  measurement  receiver. 


By  trigonometric  identities,  the  outputs  of  the  sum  and  difference 
operations  in  the  figure  are  found  to  be 
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s2(t)  = 

|^2  A  ( t)  cosy— yj 

COS 

0)Qt  +  4>(t) - 2~ 
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v' 

s3(t)  = 

2A( t)  siny— 2~ y 

sin 

V  +  *<t)  -  — 

for  most  of  the  time.  The  envelope  detectors  produce  the  magnitudes  of 
the  first  factors  in  equations  (104)  and  (105).  These  signals  pass 
through  logarithmic  amplifiers,  and  the  difference  is  taken.  Thus,  if 
Uq6  <  ir,  the  processor  input  is  proportional  to 

s4(t)  =  log  tan  u)Q6  .  (106) 

Since  the  modulation  effects  have  been  removed  and  6  is  known,  the 
processor  can  calculate  an  estimate  of  fg  =  a>g/2Tr.  The  tangent 
function  has  an  unambiguous  inverse  only  over  a  range  of  ir  radians. 
Thus ,  for  unambiguous  operation  over  the  frequency  range  of  W  hertz ,  we 
must  have  6  £  1/2W.  A  major  problem  with  the  IBM  receiver  is  that  the 
resolution  is  inversely  proportional  to  6.  Since  6  must  be  sufficiently 
small  so  that  equation  (103)  is  valid  for  most  of  the  time,  the  resolu¬ 
tion  may  be  inadequate.  Another  problem  i?  that  the  IFM  receiver  cannot 
handle  two  or  more  simultaneously  intercepted  signals  of  comparable  mag¬ 
nitudes  . 


3.5  Scanning  Superheterodyne  Raceiver 

Figure  11  shows  a  block  diagram  of  a  realization  of  a  scanning 
superheterodyne  receiver  for  frequency  estimation.  To  explain  the  oper¬ 
ation,  we  consider  the  system  response  to  one  scan  of  the  local  oscil¬ 
lator  and  an  input  that  has  constant  amplitude,  frequency,  and  phase 
over  the  scan  period,  T.  The  input  is  represented  by 

s(t)  =  A  cos  +  o)»0£t£T,  (107) 

where  w0  is  the  carrier  frequency  and  9  is  the  phase  angle  at  t  =  0, 

The  scanning  waveform,  which 
is  proportional  to 

0  <  t  <  T  ,  ( 108) 


which  defines  the  beginning  of  the  scan, 
is  the  output  of  a  swept  local  oscillator, 

y(t)  =  cos  KOgt  -  TfMt2^  , 
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Figure  11.  Scanning  superheterodyne  receiver. 


where  uis  is  the  frequency  at  t  =  0  and  u  is  the  scan  rate,  which  is  the 
rate  of  frequency  change.  The  output  of  the  mixer,  s.j(t)  =  s(t)y(t)  , 
passes  through  a  bandpass  filter  with  impulse  response  h(t)  and  band¬ 
width  2B.  Ignoring  a  high-frequency  term  that  is  suppressed  by  the 
bandpass  filter,  we  have 

S1  (t)  =  2  A  C°S  (W1fc  +  7FPt2  +  6)  '  0  -  t  -  T  '  (109) 

where  w  =  u  -  w  . 

1  0  s 

The  symmetrical  bandpass  filter  has  transfer  function  H('o)  that 
can  be  written  as 

H(w)  =  H1  (u  +  ,  (110) 

where  H^oj)  is  the  transfer  function  of  a  lowpass  filter  and  <*>c  is  the 
center  frequency  of  the  bandpass  filter.  The  first  term  on  the  right 
side  has  significant  values  only  for  positive  frequencies,  while  the 
second  term  has  significant  values  only  for  negative  frequencies. 
If  hj(t)  is  the  impulse  response  of  the  lowpass  filter,  then 

h(t)  =  2h„(t)  cos  w  t  .  (Ill) 

i  c 

The  output  of  the  bandpass  filter  is 

00 

s2(t)  =  /  s^-HhU  -  T)  dT  .  (112) 
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By  using  equations  (109)  and  (111)  and  the  pertinent  trigonometric 
relations,  equation  (112)  becomes 


s2(t)  =  J  /  h^t  -  T)  cos  |(wi  -  wc)  T  +  ^ux2  +  0  +  0>ct]  dT 

L  C  (113) 

T 

+  j/q  yt  -  X)  COS  x  +  nyx2  +  0  -  v]  dT  . 

It  is  assumed  that  Hj(u>)  has  a  sufficiently  narrow  bandwidth 
that  the  second  integral  on  the  right  side  of  this  equation  is  negli¬ 
gible.  The  time-frequency  diagram  of  figure  12,  in  which 

f  =  o)/2t r,  illustrates  the  effect  of  the  filter.  Ihe  filter  out¬ 
put,  S2(t),  is  significant  only  over  a  portion  of  the  scan  period. 

Thus,  we  can  extend  the  limits  of  the  first  integral  to  ±«°  with 

negligible  error  if  f.,  +  yT  >  fc  +  B/2  and  f1  <  f  -  B/2.  Under  these 

assumptions, 


V*’  ~iL  h,<t  '  Tl 


COS 


(V  + 


iryx2  + 


dx 


(114) 


where  U2  =  W0  -  Ws  -  ^c  an<^  =  ®  +  wct.  ^  ^urt'ler  simplify  equation 
(114),  we  assume  a  Gaussian  bandpass  filter;  that  is. 


x - H 


XIME  > 

Figure  12.  Time-frequency  diagram  for  scanning 
superheterodyne  receiver. 


(115) 
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where 


4>(  t) 


4ira4lJ(t  -  <S)2  +  4a4w2(t  -  6)  -  ttuw2 
4  ( a4  +  ir2P2) 


s3(t) 


+ 


-1/4 

exp 


2  Try6 


4^a4  +  it2!*2) 


(121) 


As  indicated  in  figure  11,  S2(t)  is  applied  to  an  envelope  detector 
that  extracts  the  envelope,  83(f) ,  from  the  input.  The  peak  value 
of  S3(t)  is  attained  when  t  =  6  -  ( o>2/2ttu )  =  6  -  (fQ  -  fg  -  fc)/V 
Thus,  the  input  frequency  f^  can  easily  be  estimated  from  the  time 
location  of  the  peak  value.  The  normalized  peak  value,  a,  which  is 
defined  as  the  peak  value  relative  to  A/2,  the  peak  value  for 
small  y ,  is 


a  = 


(122) 


The  half-power  points  of  s3(t)  are  determined  by  setting  the 
exponential  factor  in  equation  (121)  equal  to  1//2.  The  pulse  duration 
of  s3(t)  between  half-power  points  is  determined  to  be 


aI2._l_n_  .2  li^2 

irp 


(123) 


The  frequency  resolution  in  hertz.  A,  as  determined  by  the  processor  of 
figure  11,  is  approximately  equal  to  UTp,  the  frequency  range  scanned 
during  pulse  duration  Tp.  Thus,  the  resolution  is 


A  _  a(2  In  2)  1/2 

if 


(124) 


From  equation  (115),  the  3-dB  power  spectrum  bandwidth  in  hertz 
is  related  to  parameter  a  by 


_  (2  In  2)  1/2 

it  a  * 


(125) 
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In  terms  of  B,  we  can  write 


-lA 


a  =  1+0.195-^ 

V  b4> 


(126) 


A  =  B  [l  +  0.195  ) 

\  B  7 


1/2 


(127) 


These  equations  were  originally  derived  for  the  theory  of  spectrum 
analyzers . 9 

If  the  scan  rate,  P,  is  high. 


a  =  1.5  —  ,  A  =  0.44  -  ,  P  »  B2  ,  (128) 

/n  B 


which  clearly  show  the  effects  of  increasing  P. 

Using  elementary  calculus,  we  determine  the  optimal  filter 
bandwidth,  BQ ,  to  minimize  A.  Substituting  this  value  of  BQ  into 
equation  (127),  we  obtain  AQ,  the  minimum  resolution  as  a  function  of 
P.  The  results  are 


B  =  0.664/p  ,  A  =  /2  B  .  (129) 

0  0  0 


The  corresponding  normalized  peak  value  is 


(130) 


which  is  no  longer  a  function  of  P.  if  the  optimal  bandwidth  is  used, 
these  equations  indicate  that  the  achievable  resolution  becomes  worse  as 
P  increases,  but  the  peak  value  does  not  change. 


9M.  Engelson  and  F.  Telewski,  Spectrum  Analyzer  Theory  and 
Applications ,  Ar tech  House,  Inc.,  Dedham,  HA  (1974). 


a0  =0.84  , 
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Improved  resolution  at  a  high  scan  rate  can  be  achieved  by  a 
microscan  receiver ,  which  uses  the  compression  of  pulses*  The  microscan 
receiver  (fig.  13a)  includes  a  dispersive  filter ,  whereas  the  scanning 
superheterodyne  receiver  includes  a  bandpass  filter* 
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Figure  13.  Microscan  receivers  for  (a)  frequency  estimation  or 
magnitude  of  Fourier  transform  and  (b)  real  part  of 
Fourier  transform. 

The  impulse  response  of  the  dispersive  filter  is  modeled  as 


where  is  the  duration  of  the  impulse  response  and  the  amplitude  has 
been  normalized  to  unity.  The  dispersive  filter  can  be  realized  by  a 
surface  acoustic  wave  device.  10  Substituting  this  equation  and  equation 
(109)  into  equation  (112)  yields  the  dispersive  filter  response  to  one 


1UA.  A.  Oliner,  ed . ,  Acoustic  Surface  Waves 
York,  Inc.,  New  York  (1978). 


scan  of  the  local  oscillator  and  an  input  represented  by  equation  (107) 
over  the  scan  period.  Initially,  we  set  Td  =  T.  Using  2  cos  u  cos  v  = 
cos  (u  -  v)  +  cos  (u  +  v),  we  may  write  the  result  as  the  sim  of  the  two 
integrals.  We  assume  that  lw0  "  ws  ~  “c1  «  I  %  -  w  +  w  |  so  that,  for 
most  values  of  t,  we  may  neglect  one  of  the  integrals.  We  are  left  with 
the  approximation 


s  (t)  =  j  /  2  cos 
2  4 

T 

1 


(y  +  e2^  dT,  0  <  t  <  2T  , 


where 


S  (t)  =  0  ,  t  <  0  or  t  >  2T  , 


=  w  -  a)  -  w  +  2TTUt  . 
3  0  s  c 


(132) 


0  =  0  +  u)  t  -  irut2 

2  c 


(133) 


T1  =  max  (t  -  T,  0)  , 


T2  =  min  ( t,  T)  . 


If  u>3  *  o,  equation  (132)  yields 
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(134) 


The  final  form  of  equation  (134)  is  valid  even  if  w3  =  0.  For  practical 
values  of  the  parameters,  the  cosine  factor  varies  much  more  rapidly 
with  time  than  other  factors.  Consequently,  using  equations  (133),  the 
output  of  the  envelope  detector  is 


s  (t)  =7^  sine 
3  4 


,. .  A(2T  -  t) 

s  (t)  =  - - -  sine 

3  4 


(f0  -  V  fc)  '  *  ut2]  ■  0  <  1  <  *  ■ 
[('»  -  -  fc  *  “')  (2T  -  ')] 


(135) 

,  T  <  t  <  2T  , 
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where  we  have  set  u)  =  2irf 


The  peak  value  of  s3(t)  is  AT/4. 


We  define  parameter  e  as 


e  =  f  +  f  -  f.  -  pT  . 
s  c  0 


(136) 


If  the  frequency  of  the  intercepted  signal  is  such  that  e  =  0,  then  the 
peak  value  of  s3(t)  occurs  at  t  =  T;  that  is, 

s  (T)  =  7^  ,  e  =  0  .  (137) 

3  4 

We  also  have 


These  equations  indicate  that  1/pT  is  an  approximate  measure  of  the 
width  of  the  compressed  output  pulse.  Satisfying  the  inequality  ensures 
that  the  response  of  the  microscan  receiver  due  to  one  scan  does  not 
interact  significantly  with  the  response  due  to  the  next  scan. 


If  the  frequency  of  the  intercepted  signal  shifts  slightly, 
in  e  t  0.  A  small  shift  yields 


l£l  «  T  , 
V 


(139) 


where  the  right  side  of  the  equation  is  the  peak  value  of  s3(t).  Thus, 
the  peak  value  occurs  at  time  t  =  T  +  e/p.  Using  equation  (136),  we  can 
estimate  the  input  frequency,  fQ  =  u^/2i;,  from  the  time  location  of  the 
peak  value. 


We  define  the  resolution  of  the  microscan  receiver.  A,  as  the 
value  of  the  frequency  shift  in  hertz,  e,  that  produces  a  shift  in  the 
time  of  the  peak  output  equal  to  the  width  of  the  output  when  e  =  0. 
Thus, 


(140) 


where  W  =  pT  is  the  total  bandwidth  scanned 


The  ratio  of  the  input  pulse  width  to  the  compressed  output 
pulse  width  is  TW,  which  is  called  the  compression  ratio  of  the  micro¬ 
scan  receiver.  When  this  ratio  is  large,  a  comparison  of  equations 
(140)  and  (129)  indicates  that  there  is  a  substantial  improvement  in 
resolution  of  the  microscan  receiver  over  the  scanning  superheterodyne 
receiver.  A  comparison  of  the  receiver  outputs  for  typical  inputs  of 
the  form  of  equation  (107)  is  depicted  in  figure  14. 


Figure  14.  Response  of  scanning  superheterodyne  and 

irucroscan  receivers  to  simultaneous  signals. 

Due  to  the  presence  of  the  sine  function,  equation  (135) 
exhibits  smaller  subsidiary  peaks  in  addition  to  the  main  peak.  In  a 
multiple-target  environment,  these  side  lobes  can  mask  adjacent  main 
peaks  due  to  other  targets,  thus  limiting  the  frequency  resolution  and 
dynamic  range  of  the  receiver.  Consequently,  it  is  sometimes  desirable 
to  follow  the  dispersive  filter  with  a  shaping  filter  for  side-lobe 
reduction  or  to  modify  the  impulse  response  of  the  dispersive  filter 
itself . 


For  a  large  compression  ratio,  the  transfer  function  associated 
with  equation  (131)  has  a  flat,  nearly  rectangular  amplitude  response 
and  a  quadratic  phase  response11  over  bandwidth  W. 

The  preceding  analysis  is  valid  if  the  modulation  period  of  the 
input  is  large  compared  with  the  scan  period.  If  a  more  rapidly  modu¬ 
lated  input  is  present,  the  microscan  receiver  can  be  designed  to 

llH.  J.  Blinchikof  f  and  A.  I.  Zverev,  Filtering  in  the  Time  and 
Frequency  Domains,  John  Wiley  and  Sons,  Inc.,  New  York  (1976). 
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produce  an  output  that  is  an  approximation  of  the  Fourier  transform  of 
the  modulation.  In  this  case,  the  system  is  often  called  a  chirp  trans¬ 
form  processor.  As  previously,  we  consider  the  system  response  to  one 
scan  of  the  local  oscillator.  Over  the  scan  period,  the  input  is 
assumed  to  have  the  form 


i  ( t )  =  s  ( t )  cos  (w  t  +  ,  0<t< 

0  Vo  )  ~  ~ 


(141) 


This  input  is  mixed  with  the  scanning  waveform  of  equation  (108)  to 
produce 


;  (t)  =  —  s  (t)  cos  (to  t  +  init2  +  s\  /  0£t£T. 

1  2  0  V  1  / 


(142) 


The  impulse  response  of  the  dispersive  filter  is  given  by  equation 
(131),  where  we  assume  that  T^  >_  T.  we  calculate  the  system  output  in 
the  time  interval  T  _<  t  <_  T^.  Substituting  equations  (142)  and  (131) 
into  equation  (112),  using  trigonometry,  dropping  a  negligible  integral, 
and  substituting  a  complex  exponential,  we  obtain 


s 2  ( t )  =  Re  <4  exp 


=  Re  exp 


(“cfc  -  sQ(T)q(T)exp|-jT  ^o)c  -  co ^  -  2Tryt)J  dfj 

£  (“ct  -  *pt2+ejj  SQ  (u)c  -  o)1  -  2irUt)| 


=  J  |s  (u>c  -  0)  -  2Trut^  |  cos  [V  -  irut2+0  +  <t> 


0)^-0)  -  2TT)jtj 


T  <  t  <  T .  , 
—  —  a 


(143) 


where  q(t)  is  defined  in  equation  (85),  Sq(o))  is  the  Fourier  transform 
of  Sg(T)q(T),  and  4>(<d)  is  the  phase  angle  of  the  Fourier  transform.  The 
output  of  the  envelope  detector  is  proportional  to 


s  (t)  =  |S  /co  + 
3  0  v  s 


0)  -  0) 
c  0 


2irjjt)  |  ,  T  <  t  <  Td  . 


(144) 


Thus,  the  magnitude  of  the  Fourier  transform  of  the  input 
modulation  has  been  produced  as  a  time  signal. 


45 


If  s2(t)  is  multiplied  by  h(t),  then  equations  (131)  and  (143) 
indicate  that,  after  elimination  of  a  double  frequency  term  by  filter¬ 
ing,  the  phase-shifted  real  part  of  Sq(u)  is  produced  as  a  time  sig¬ 
nal.  As  shown  in  figure  13(b)  the  waveform  h(t),  T  _<  t  <_  Td,  can  be 
produced  by  applying  an  impulse  at  time  T  to  a  dispersive  filter  with 
impulse  response  h.,(t)  =  h(t  +  T) ,  0  _<  t  <_  Td  -  T.  The  phase-shifted 
imaginary  part  of  Sg  ( u>)  can  be  produced  as  a  time  signal  by  multiplying 
s2(t)  by  sin  (wct  -  irpt2). 

Let  um  denote  the  maximum  value  of  u>  for  which  |  Sg  ( cu )  |  has  a 
significant  value.  If 


U)  Cl)  +  Cl) 

_  ,  m  _  s  c 

T  +  —  <  - r — 

2v\i  —  2ttv 


lTd  ' 


(145) 


then  during  the  time  interval  T  <_  t  <_  Tj,  s3(t.)  exhibits  all  the  values 
of  |Sg(u>)|  for  0  <_  u)  <_  wm.  The  range  of  possible  values  of  fg  =  u)q/2tt 
for  which  equation  ( 145)  can  be  satisfied  is 


W„  =  pT,  -  f  -  pT 
R  d  m 


T  >  T  + 
a  — 


(146) 


where  fffl  =  tom/2Tr.  If  we  wish  to  avoid  the  interference  of  the  Fourier 
transform  generated  by  a  scan,  which  occurs  during  T  <_  t  £  T  ,  with  the 
Fourier  transform  generated  by  the  next  scan,  then  we  set  T^  £  2T.  In 
terms  of  the  total  bandwidth  scanned,  W  =  pT,  we  have 


0<W<w(l--p),  T  +  ——  <  T  <  2T  .  (  147) 

—  R—  \  W  /  p  —  d  — 

If  fm  «  W  and  Td  =  2T,  then  WR  =  W. 

The  inequality  for  Tj  in  equation  (147)  is  necessary  for 
satisfactory  spectral  analysis.  However,  for  frequency  estimation 
alone,  setting  Tj  =  T  not  only  is  adequate,  but  also  minimizes  the 
interference  between  scan  responses. 

The  chirp  transform  processor  can  be  used  as  the  basic  building 
block  of  an  analog  version  of  the  cross  correlator  of  figure  5. 

In  addition  to  frequency  estimation,  the  microscan  receiver  can 
be  used  for  detection.  Alternatively,  if  interception  is  verified  by  a 
parallel  system  designed  expressly  for  that  purpose,  frequency  can  be 
estimated  by  the  microscan  receiver. 
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If  the  scanning  period  is  less  than  the  period  between 
frequency  changes  of  the  intercepted  signal,  the  microscan  and  scanning 
superheterodyne  receivers  can  estimate  each  frequency  of  a 
frequency-hopping  or  MFSK  signal.  If  not,  some  frequencies  may  be 
missed,  and  the  estimation  accuracy  is  degraded.  The  inherent  linearity 
of  these  two  receivers  makes  them  potentially  effective  in  analyzing 
many  simultaneous  intercepted  signals. 


4.  DIRECTION  FINDING 

Signals  must  be  detected,  and  sometimes  the  frequency  must  be 
estimated,  if  direction  is  to  be  found.  Conversely,  direction  finding 
provides  signal  sorting,  which  restricts  the  number  of  signals  that  the 
detection  and  frequency  estimation  systems  must  process  simultaneously. 

For  simplicity,  the  estimation  of  a  single  bearing  angle  is 
considered.  In  a  ground-based  interception  system,  an  azimuth  angle  may 
be  all  that  is  needed.  However,  airborne  systems  may  require  estimates 
of  both  the  azimuth  and  the  elevation  angles  to  the  intercepted 
transmitter . 

4.1  Energy  Comparison  Systems 

Energy  comparison  systems  are  analogous  to  the  amplitude 
comparison  systems  used  in  radar ^ ,  but  when  communications  with  unknown 
parameters  are  to  be  intercepted,  it  is  logical  to  base  comparisons  upon 
the  energy  rather  than  the  amplitude.  The  stationary  multibeam  system 
for  direction  finding  is  illustrated  in  figure  15.  First,  the  largest 
of  the  receiver  outputs  is  selected.  Next,  the  larger  of  the  two 
receiver  outputs  corresponding  to  beams  adjacent  to  the  beam  that  pro¬ 
duced  the  largest  output  is  selected.  The  two  selected  outputs  are 
denoted  by  and  L2«  The  processor  compares  the  outputs  to  a  threshold 
for  detection.  The  angle  of  arrival  is  estimated  from  the  logarithm  of 
the  ratio  of  L1  to  L2>  The  radiation  patterns  of  the  adjacent  beams  are 
illustrated  in  figure  16,  where  <f>  represents  the  angle  of  arrival  of  an 
intercepted  signal,  and  beam  pattern  Fi(6)  produces  L^.  The  origin  of 
the  coordinate  system  is  defined  so  that  +\|»  and  -\J>  indicate  the  peak 
responses  of  the  two  beams,  respectively.  Suppose  that  the  beam 
patterns  are  approximately  Gaussian;  that  is,  they  are  described  by 


F  (0)  =  K1  exp 


(  148) 


12P.  K.  Barton, 
(1976). 


Radar  System  Analysis ,  Artech  House, 
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f2(B)  =  k2 


exp  - 


(  149) 


w^ere  and  K2  are  constants  independent  of  0,  and  b  is  a  measure  of 
the  beam  width.  If  the  receivers  contain  radiometers,  then  the  Li  are 
proportional  to  the  squares  of  the  F^.  Thus,  in  the  absence  of  noise, 
the  processor  input  is 


h2(*]  K2  b2 


(150) 


The  angle  of  arrival  can  be  determined  by  inverting  this  equation.  In 
the  presence  of  noise,  the  same  inverse  provides  an  estimate,  of 

the  actual  angle  <f>.  We  have 


(151) 


where  Z  is  now  a  random  variable.  Since  i  is  a  linear  function  of  Z, 
the  required  processing  is  quite  simple.  We  give  an  error  analysis 
assuming  that.  and  are  correctly  chosen. 
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Figure  15.  Stationary 
multibeam  system. 


Figure  16.  Adjacent  antenna 
radiation  patterns. 


Because  L1  and  L2  are  the  outputs  of  radiometers,  they  are 
statistically  independent,  noncentral  y2  random  variables  with  y  degrees 
of  freedom  and  noncentral  parameters  X1  and  X2,  respectively.  From 
elementary  probability  theory,  the  probability  density  function  of  Y  = 
L^/L 2,  the  quotient  of  two  nonnegative  random  variables,  is 


p(y)=  /  xp  (yx)p  (x)  dx  ,  y>0 

2  0  11  12 


p  (y)  =  0  ,  y  <  0 

2 


(152) 


where  p^fx)  and  p12(x)  are  the  density  functions  for  L1  and  L2, 
respectively.  The  density  function  p^fx)  is  given  by  equation  (30) 
with  X^  substituted  for  X.  The  noncentral  parameters  are  given  by 


2E 


1 


x,  -  ^  exp 


_  2(d)  -  »)2 
b2 


2E2 


2(d>  +  d))2 
b2 


(153) 


where  E-  is  th.i  energy  received  when  the  intercepted  signal  enters  the 
center  of  the  beam  associated  with  . 


From  equations  (151)  to  (153)  and  the  fact  that  Z  =  In  Y,  we 
can  determine  the  root-mean-square  error  of  $,  which  we  denote  by  ER. 
By  definition  and  a  straightforward  expansion. 


E2  =  E  [($  -  d>)2]  =  VAR(<ji)  +  B2  , 


where  VAR($)  is  the  variance  of  $  and  Bx  is  the  bias. 


(154) 


(155) 


To  evaluate  ER,  we  need  the  first  two  moments  of  Z.  For  simplicity,  we 
assume  that  K1  =  K2  and  E1  =  E2<  When  d>  =  0,  equations  (153), give  X1  = 
X2  so  that  L1  and  L2  are  identically  distributed  random  variables. 
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Since  Z  =  In  -  In  L2'  it  follows  that 

E[Z]  =  0  ,  (ft  =  0  .  (156) 


At  other  values  of  closed-form  expressions  for  the  mean  and  other 
moments  are  difficult  to  obtain.  In  general,  the  distribution  of  L1  at 
<)>  is  identical  to  the  distribution  of  L2  at  -$.  Thus,  ER  is  symmetric 
about  <j>  =  0 . 

To  obtain  an  approximate  expression  for  E[Z],  we  write  the 
integral  for  E[ln  Y]  as  a  double  integral  and  change  coordinates: 


00  00 


E[Z]  =  /  /  x  In  y  p  (yx)p  (x)  dx  dy 

00  11  12 


00  00 

=  /  /  ln/-\p  ( v ) p  (x)  dx  dv 

0  0  W  11  12 


(157) 


The  logarithm  is  approximated  by  the  first  six  terms  of  its 
two-dimensional  Taylor  series  expansion  about  the  point  v  =  m^ ,  x  =  m2, 
where  m1  and  m2  are  the  mean  values  of  L1  and  L 2,  respectively.  Thus, 


In 


e> 


v  -  x  - 

in  .  m  _ 


(v  -  mi  y 


2m2 


+  (X  ~  m2)2  _  (v  ~  ml)(X  ~  m2) 


(158) 


2m| 


2m^  m^ 


This  approximation  is  accurate  over  some  range  of  v  and  x  about  v  =  m^, 
x  =  m2.  If  p-j-jfv)  and  p12(x)  are  negligible  outside  this  range,  then 
the  substitution  of  equation  ( 158)  into  equation  ( 157)  yields  an  accu¬ 
rate  approximation  of  E[Z).  The  ranges  of  significant  values  of  p^fv) 
and  p^2(x)  are  approximately  limited  by  |v  -  m^|  <  3o^  and  |x  -  m2l  < 
3o2,  where  and  a2  are  the  standard  deviations  of  and  L2,  respec¬ 
tively.  Thus,  sufficient  conditions  for  the  valid  use  of  equation  (158) 
are  <<  m^  and  ct2  <<  m2. 
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Making  the  substitution  and  using  the  properties  of  density 
functions,  we  obtain 


/mA  °i  °2 

E  [Z]  =  lnl — ]-  -  +  -  »  o,  <<  m,  ,  a„  <<  m„  . 

\2/  2m2  2m2  1  1  2  2 


(159) 


From  equations  (26)  and  (27),  we  have 


“1  -  xi  + 

2y  ,  i  =  1,  2  , 

(160) 
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*3 

(161) 

An  approximate  expression  for  E[Z<-]  is  obtained  in  an  analogous 
manner.  Combining  this  expression  with  equation  (159)  and  dropping 
terms  higher  than  second  order  in  and  02/1*12*  we  obtain 


VAR(Z) 


°1  «  ml 


d„  << 

2  2 


(162) 


By  using  K1  =  K2  and  equations  (151),  (154),  (159),  and  (162),  an 

equation  for  ER  can  be  derived, 


a 

_ 2 _ 8ij^ 

2m2  b2 


°1  «  ",  / 


<<  m2  * 


(163) 


Figures  17  and  18  show  representative  plots  of  ER  versus  4>. 
Figure  17  illustrates  the  effect  of  the  normalized  beam  width,  b/i|>,  on 
ER.  For  E/Ng  =  lO*4,  y=  10^,  and  p  =  p  the  normalized  beam  width  that 

minimizes  is  approximately  4.5.  However,  the  curves  show  that  this 

normalized  beam  width  is  not  optimal  for  other  values  of  <j>.  Figure  18 

illustrates  the  effect  of  increasing  the  observation  interval,  T.  Since 

an  increase  in  T  causes  proportionate  increases  in  E  and  Y,  we  set  E/Ng 
=  10y  for  each  value  of  Y«  The  optimal  normalized  beam  width  for  each 
value  of  Y  and  <j>  =  i()  is  chosen.  The  curves  show  a  steady  decrease  in  ER 
as  Y  increases. 
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ARRIVAL  ANGLE.  o/y, 


ARRIVAL  ANGLE.  o/y, 


Root-mean- square  error  versus  arrival  angle 
for  optimal  beam  widths  and  different  values 
of  y. 


Suppose  many  frequency-hopping  signals  are  simultaneously 
intercepted  by  a  channelized  receiver.  If  it  is  impossible  to  correlate 
successive  hopping  frequencies,  the  direction  of  a  signal  must  be  esti¬ 
mated  on  the  basis  of  the  energy  received  in  a  single  channel  during  a 
single  hopping  period.  With  this  interpretation,  y  is  proportional  to 
the  hopping  period,  and  figure  18  indicates  the  loss  in  directional 
accuracy  resulting  from  an  increase  in  hopping  rate. 


Energy  comparison  with  a  rotating  beam  is  illustrated  in  fig¬ 
ure  19.  A  rotating  dish  or  an  electronically  scanned  phased  array  is 
used  with  an  omnidirectional  antenna  a  small  distance  apart.  The 
receivers  are  radiometers.  Detection  of  a  signal  is  verified  and  an 
estimate  of  direction  is  obtained  by  measuring  the  ratio  of  L1  to  L2  and 
comparing  the  ratio  to  a  threshold.  Because  a  ratio  is  used,  the 
effects  of  amplitude  modulation  and  rotation  can  be  largely  eliminated 
by  the  processor.  Signals  entering  through  the  side  lobes  and  the  back 
lobes  of  the  rotating  beam  are  inhibited.  Alternatively,  if  two  simul¬ 
taneously  rotating  beams  are  offset  relative  to  each  other,  sum  and 
difference  outputs  can  be  produced  (fig.  20),  as  in  monopulse  radar 
systems.  If  the  intercepted  signal  duration  is  sufficient,  the  accuracy 
of  the  direction  estimates  of  the  systems  of  figures  19  and  20  can  be 
improved  by  continued  processing.  After  the  target  is  detected,  the 
rotation  slows  or  stops.  Processing  over  an  increased  observation 
interval  allows  further  adjustment  of  antenna  position  until  the  inter¬ 
cepted  signal's  direction  of  arrival  is  fixed  near  the  center  of  the 
rotating  beam  or  beams. 


Figure  19.  Rotating  beam  system. 


Figure  20.  Rotating  monopulse  system. 


The  performance  analysis  of  the  system  of  figure  19  is  similar 
to  the  analysis  of  the  stationary  multibeam  system.  We  give  an  analysis 
of  the  monopulse  system. 

Assuming  that  the  receivers  of  figure  20  contain  radiometers, 
is  proportional  to  the  square  of  the  pattern  ( <J>) .  It  follows  that, 
in  the  absence  of  noise. 


L1  "  L2  F1(<(,)  "  F2(4>) 

L  +  I.  _  ^  J  ‘ 

1  2  Fi(*)  +  p2(*) 


Y 


(164) 


At  <(>  =  0,  we  have  m.,  =  m2  =  m  and  o1  =  cs2  =0  so  that 


E  -  ^b-2q 
R  8>Jm  ' 


<(>  =  0  ,  a  <<  in  . 


(169) 


Appropriate  substitution  into  this  equation  indicates  that  we  can 
achieve  ER  <<  b  if  the  observation  interval,  T,  is  sufficiently  large. 
Thus,  the  potential  direction-finding  accuracy  is  much  better  than  a 
beam  width. 

To  compare  this  result  with  the  corresponding  result  for  the 
stationary  multibeam  array,  we  set  m.j  =  m2  =  m,  =  0 2  =  and  41  =  0 

in  equation  (163).  The  result  is  equation  (169).  Thus,  if  <j>  =  0,  the 
stationary  multibeam  array  performs  as  well  as  the  monopulse  system. 
However,  if  *  0,  ER  increases  in  the  multibeam  case.  In  the  monopulse 
case,  if  the  intercepted  signal  duration  is  sufficient,  the  antennas  can 
adjust  their  positions  until  <f>  -  0  so  that  the  ER  after  adjustment  is 
given  by  equation  (169).  The  main  disadvantage  with  the  monopulse  or 
other  rotating  beam  system  is  the  narrow  instantaneous  field  of  view, 
which  may  cause  a  signal  to  be  missed  or  may  decrease  the  possible 
observation  time.  On  the  other  hand,  when  many  hostile  communications 
are  present,  the  narrow  field  of  view  provides  a  valuable  signal-sorting 
capability. 

4. 2  Interferometer 

Another  direction-finding  system  is  the  interferometer ,  which 
also  forms  the  heart  of  phase- comparison  monopulse  radar  systems.1^  An 
interferometer  consists  of  two  or  more  antennas  or  groups  of  elements  of 
a  phased  array  that  use  phase  or  arrival-time  information  to  estimate 
direction  (fig.  21).  The  antennas  may  be  rotating  and  directional. 

We  first  consider  the  interferometer  of  figure  21(a).  Suppose 
a  plane  wave  arrives  at  angle  <fi,  where  I  <f  I  <_  t/2.  If  two  antennas  are 
separated  by  distance  d,  then  since  phase  angles  are  modulo  2n  numbers, 
the  phase  difference  between  the  antenna  outputs  is 


e  =  2”d  ?in-i  _  2  nn  ,  1 9 1  <  ir  ,  |<fr|  <  J  ,  (170) 


120.  K.  Barton,  Radar  System  Analysis ,  Ar tech  House,  Inc.,  Dedham,  MA 
(1976). 


•  * 


where  Ag  is  the  signal  wavelength  and  n  is  an  integer  that  ensures 
satisfaction  of  the  first  inequality.  If  ci  <_  As/2,  then  n  =  0.  if 
d  >  Ag/2,  then  n  varies  with  taking  negative  and  positive  values  and 
the  value  zero , 

PLANE  WAVE 


(») 


Figure  21.  Interferometers  using  (a)  phase  information  and 
(b)  arrival -time  information. 


One  antenna  output  provides  a  reference.  Each  of  the  other 
antenna  outputs  is  applied  to  a  separate  device  that  estimates  its  phase 
relative  to  the  reference.  This  device  nay  be  similar  to  the  part  of 
the  IFM  receiver  in  figure  10  that  is  fed  by^s(t)  and  s^(t).  The  esti¬ 
mates  of  the  relative  phases  are  denoted  by  0^. 
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„  Estimates  of  the  angle  of  arrival  are  calculated  from  the 

0^  by  inverting  equation  (170).  If  there  are  m  antennas,  the  (m  -  1) 
estimates  are  given  by 


(171) 


I't’jJ  <  \  ,  ISjJ  <  ",  i 


2 ,  .  .  . ,  m- 1 , 


where  A„  is  an  estimate  of  . 

s  s 

Unless  d^  Ag/2,  tbe  appropriate  value  of  n^  is  unknown.  In  fact, 
many  different  values  of  n^  may  satisfy  the  first  inequality.  On  the 
other  hand,  the  accuracy  of  the  estimate  of  4>  increases  with  d^.  To 
obtain  an  estimate  of  <)>  that  is  both  accurate  and  unambiguous,  three  or 
more  antennas  are  used.  The  first  antenna  is  a  reference.  The  second 
antenna  is  separated  by  d^  _£  Ag/2  from  the  reference.  The  other 
antennas  are  separated  by  increasingly  greater  distances  from  the 
reference.  The  first  antenna  pair  provides  an  unambiguous  estimate  of 
the  angle  of  arrival,  <f>.  The  first  and  third  antennas  provide  a  more 
accurate  estimate  of  tj>  but  an  ambiguous  one.  However,  this  ambiguity  is 
resolved  by  the  first  estimate.  Subsequent  antennas  allow  increasingly 
accurate  estimates  of  <J>,  provided  that  the  ambiguities  can  be  resolved 
by  the  less  accurate  estimates.  In  designing  the  interferometer, 
distance  d1  should  correspond  to  the  shortest  wavelength  to  be 
intercepted . 

Assuming  that  n^  is  correctly  determined,  the  probability 
density  function  and  moments  of  ^  can  be  calculated  by  standard 
methods  once  the  joint  probability  density  function  of  0^  and  As  is 
specified.  If  there  is  a  negligible  error  in  As,  only  the  probability 
density  function  of  0^  is  required.  Both  of  these  density  functions 
depend  on  the  details  of  the  receiver  implementation. 


An  alternative  method,  depicted  in  figure  21(b),  which 
eliminates  the  need  to  estimate  As  and  the  need  to  resolve  ambiguities, 
is  based  on  the  direct  measurement  of  the  relative  arrival  time  of  a 
plane  wave  at  two  antennas.  This  relative  arrival  time  can  be  estimated 
by  using  the  vector  output,  ,  of  the  cross  correlator  in  figure  5,  or 
by  using  two  arrival-time  estimators. 


Assuming  that  a  plane  wave  is  received,  the  relative  arrival 
time  is  given  by 


where  v@  is  the  speed  of  electromagnetic  waves.  Consequently,  the  angle 
of  arrival  estimator  is 


(173) 


where  Tr  is  the  relative  arrival-time  estimator.  Retaining  the  first 
three  terms  of  a  Taylor  series  expansion  of  equation  (173)  about  the 
point  Tr  and  using  equation  (172),  we  obtain 


je_ 


d  cos  ij> 


v|  sin  <)> 
2d^  cos  ^  <(> 


l<H  <  it/2  ,  (174) 


If  Tr  provides  an  unbiased  estimate  of  Tr,  then 

e[tJ  =  Tr  . 


It  follows  that 


(175) 


<t>  + 


v|  sin  <() 
2d^  cos^  <)> 


l<t>l  <  tt/2  . 


(176) 


Equation  (176)  shows  that  the  angle  estimator  is  biased  even  if  the 
arrival-time  estimator^  is  unbiased.  A  straightforward  calculation 
yields  the  variance  of 

The  root-mean-square  error  of  £  is  obtained  by  substituting 
equation  (174)  into  equation  (154).  By  retaining  only  the  lowest  order 
term,  the  result  is 


;R  “  d  cos  <(>  °(Tr)  '  I'M  <  2  ' 
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where  a( Tr)  is  the  standard  deviation  of  Tr<  If  <p  =  0,  ER  has  its  min¬ 
imal  value,  we  have 


This  equation  indicates  that  large  antenna  separations  are  sometimes 
necessary  if  ER  is  to  be  small.  For  example,  if  c(Tr)  =  100  ns  and  we 
desire  ER  =  0.05  radians,  then  we  need  d  =  600  m.  However,  unless  the 
antenna  separation  is  much  smaller  than  the  distance  to  the  source  of 
the  intercepted  signal,  the  plane  wave  assumption  and,  hence,  equation 
{176)  do  not  apply. 

To  determine  VAr(  Tr ) ,  we  must  specify  the  details  of  the  rela¬ 
tive  arrival-time  estimation.  For  example,  if  the  relative  arrival  time 
is  computed  by  subtracting  one  arrival  time  from  another  and  if  the 
arrival  times  are  statistically  independent  with  equal  variances, 
VAr( Tr )  is  twice  the  variance  of  the  arrival  times.  Expressions  for 
the  latter  variance  and  various  methods  for  arrival-time  estimation  can 
be  found  in  the  literature.13'14 

5.  CONCLUSIONS 

When  little  is  known  about  a  signal  to  be  intercepted,  the  signal 
energy  and  the  autocorrelation  function  present  two  natural  character¬ 
istics  upon  which  the  interception  receiver;  can  base  its  processing.  A 
radiometer,  which  measures  the  signal  energy,  and  a  cross  correlator, 
which  estimates  the  autocorrelation  function,  perform  similarly  as 
detectors.  The  choice  between  using  a  radiometer  and  using  a  cross 
correlator  is  best  made  on  the  basis  of  the  additional  hardware  needed 
for  the  receiver's  other  functions.  For  example,  if  a  Courier  analysis 
of  the  signal  and  the  direction  finding  with  an  interferometer  are 
planned,  then  detection  with  a  cross  correlator  is  preferable  since  it 
requires  little  additional  hardware. 

A  receiver  is  channelized  if  its  total  bandwidth  is  divided  into  M 
parts  by  a  filter  bank,  and  the  filter  outputs  are  processed  in  paral¬ 
lel.  Against  narrowband  communications,  the  intercepted  power  required 
by  a  channelized  receiver  is  reduced  by  a  factor  of  /it  compared  with  a 
wideband  receiver  with  the  same  total  bandwidth.  Channelized  receivers 
often  improve  performance  against  frequency-hopping  communications, 
although  a  sufficiently  high  hopping  rate  may  render  a  channelized 
design  impractical.  Pseudonoise  spread-spectrum  communications  require 
preliminary  processing  before  being  applied  to  a  channelized  receiver. 

B.  N.  Mityashev ,  The  Determination  of  the  Time  Position  of  Pulses  in 
the  Presence  of  Noise,  MacDonald ,  London  (196 5). 

J4D.  J.  Torrieri,  Adaptive  Thresholding  Systems,  TEES  Trans.  Aerosp. 
Electron.  Syst.,  1J_  (May  1977),  273. 
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Channelized  receivers,  digital  filters,  acousto-optical  devices,  and 
microscan  receivers  (chirp  transform  processors)  provide  comparable 
resolutions  as  frequency  estimation  systems.  The  systems  also  theoret¬ 
ically  can  handle  multiple  signals  that  are  simultaneously  inter¬ 
cepted.  The  IKM  and  scanning  superheterodyne  receivers  are  less 
attractive  systems  in  terms  of  potential  capabilities,  but  may  be  useful 
as  auxiliary  frequency  estimators.  If  speed  of  frequency  estimation  is 
important,  it  is  desirable  to  use  channelized  receivers,  microscan 
receivers,  or  acousto-optical  devices  with  parallel  readout  of  the 
photodetector  outputs.  All  three  systems  can  provide  detection  with 
little  additional  hardware.  In  choosing  among  these  three  systems,  cost 
and  practical  implementation  problems  are  among  the  most  important 
considerations . 

I 

If  widely  separated  antennas  can  be  deployed  and  if  the  signal  to  be 
intercepted  has  a  narrow  autocorrelation  function,  then  an  interfer¬ 
ometer  is  an  effective  system  for  direction  finding.  Under  other  cir¬ 
cumstances,  energy  comparison  systems  and  possibly  acousto-optical 
devices  are  attractive  alternatives.  Energy  comparison  systems  can  be 
either  stationary  or  rotating.  Rotating  systems  can  adjust  the  antenna 
position  to  increase  the  angle  estimation  accuracy.  When  many  hostile 
communications  are  present,  the  narrow  field  of  view  provides  a  valuable 
signal-sorting  capability.  However,  a  narrow  instantaneous  field  of 
view  may  cause  a  signal  to  be  missed  or  may  decrease  the  possible 
observation  time,  thereby  decreasing  the  angle  estimation  accuracy. 

By  what  electronic  countermeasures  can  the  communicators  thwart 
interception?  The  data  rates  and  the  transmission  powers  can  be  kept  to 
a  minimum.  Cables  and  optical  fiber  links  are  very  helpful  whenever 
feasible.  Directional  antennas  help  to  conceal  the  existence  of  commun¬ 
ications  from  the  opponent.  However,  there  are  constraints  on  the 
degree  of  directionality  that  can  be  designed  into  an  antenna  to  be  used 
in  the  battlefield.  An  important  constraint  is  the  need  to  keep  the 
antenna  small  to  hide  it  from  sight. 

Since  the  antenna  beam  angle  can  be  decreased  by  the  use  of  a 
smaller  wavelength  as  well  as  by  a  larger  antenna,  millimeter  or  even 
higher  frequencies  are  sometimes  viable  alternatives  to  radio  frequen¬ 
cies.  The  decision  to  use  smaller  wavelengths  is  tempered  by  such 
things  as  cost,  available  power,  and  propagation  properties.  The 
shorter  wavelengths  in  general  are  attenuated  more  than  the  longer 
wavelengths  and  are  more  easily  blocked  by  obstructions  in  their  path. 
Furthermore,  if  the  beam  width  is  exceedingly  narrow,  it  is  difficult  to 
keep  it  centered  on  another  station  of  a  communication  network. 

Time  hopping,  in  which  transmissions  are  increased  in  total  dur¬ 
ation,  but  contain  pseudorandom  time  gaps,  is  another  general  counter¬ 
measure.  Because  of  the  pseudorandom  gaps,  an  interception  receiver 
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must  either  process  more  noise  energy  than  it  would  otherwise  or 
decrease  its  observation  interval.  In  either  case,  performance 
degrades. 

Spread-spectrum  communications  are  inherently  more  difficult  to 
intercept  than  are  conventional  communications.  Frequency  hopping  over 
a  wide  bandwidth,  a  powerful  countermeasure,  complicates  the  design  of 
an  interception  receiver  and  degrades  its  performance.  Detection, 
frequency-estimation,  and  direction-finding  difficulties  increase  with 
the  hopping  rate. 
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GLOSSARY  OF  PRINCIPAL  SYMBOLS 

a  Parameter  proportional  to  bandwidth 

A  Signal  amplitude 

A(t)  Amplitude  modulation 

b  Measure  of  beam  width 

B  Bandwidth  of  filter 

d  Distance  between  two  antennas 

E  Signal  energy 

Er  Root-mean-square  error  of  direction  estimate 

E[x]  Expected  value  of  x 

fc  Center  frequency  of  filter  (hertz) 

fs  Frequency  of  swept  local  oscillator  (hertz) 

fp  Carrier  frequency  of  intercepted  signal  (hertz) 

F(0)  Beam  pattern 

h(t)  Impulse  response  of  filter 

H(u)  Filter  transfer  function 

In(x)  Modified  Bessel  function  of  order  n 

L  Receiver  output 

m(t)  Binary  message  sequence 

M  Number  of  channels  in  channelized  radiometer 

n(t)  White  Gaussian  noise 

N  Number  of  comparator  outputs  examinee  during  observation 

interval 

Nq  Twice  noise  power  spectral  density 

N1  Number  of  comparator  outputs  examined  during  signal  duration 

p(x)  Probability  density  function 

PD  Probability  of  detection 
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GLOSSARY  OF  PRINCIPAL  SYMBOLS  (Cont’d) 

Probability  that  particular  radiometer  exceeds  threshold 
when  signal  is  present  in  that  radiometer 

Probability  that  some  radiometer  exceeds  threshold  at  end  of 
sampling  interval  when  signal  is  present 

Probability  of  false  alarm 

Probability  that  particular  radiometer  output  exceeds 
threshold  when  no  signal  is  present 

Received  waveform,  usually  after  bandpass  filtering 

Received  power  of  intercepted  signal 

Signal  component 

Waveform  resulting  from  truncation  of  s(t) 

Fourier  transform  of  s'(t) 

Observation  interval;  period  of  single  scan 

Duration  of  impulse  response  of  dispersive  filter 

Delay  in  arrival  time  of  intercepted  signal  at  one  antenna 
output  relative  to  arrival  time  at  another  output 

Duration  of  sampling  interval  of  channelized  radiometer 

Signal  duration 

Acoustic  velocity 

Electromagnetic  velocity 

Test  statistic 

Fixed  threshold  level  of  comparator 
Variance  of  x 

Total  bandwidth  of  system  (hertz) 

Bandwidth  of  constituent  radiometers  of  channelized 
radiometer 

Normalized  peak  value  of  response  of  scanning 
superheterodyne  receiver 

Parameter  defined  by  equation  (4) 

Parameter  defined  by  equation  (44) 
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GLOSSARY  OF  PRINCIPAL  SYMBOLS  (Cont'd) 

Largest  integer  less  than  or  equal  to  TW 

Time  delay  due  to  filter  or  device 

Resolution  of  frequency  estimation  system 

Parameter  defined  by  equation  (136) 

Largest  integer  less  than  or  equal  to  TgWg 

Phase  angle;  geometrical  angle 

Noncentral  parameter  of  chi-square  distribution 
Acoustic  wavelength 
Optical  wavelength 
Signal  wavelength 

Rate  of  frequency  change  (hertz  per  second) 

Parameter  defined  by  equation  (5) 

Parameter  defined  by  equation  (45) 

Variance 

Arrival  angle  of  intercepted  signal 
Angle  modulation  function 

Chi  squared,  pertaining  to  probability  distribution 
Parameter  related  to  maximum  of  radiation  pattern 
Center  frequency  of  filter  (radians  per  second) 

Frequency  of  swept  local  oscillator  (radians  per  second) 
Carrier  frequency  of  intercepted  signal  (radians  per  second 
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